Phase synchronization on scale-free and random networks in the presence of noise 
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In this work we investigate the stability of synchronized states for the Kuramoto model on scale-free and 
random networks in the presence of white noise forcing. We show that for a fixed coupling constant, the 
robustness of the globally synchronized state against the noise is dependent on the noise intensity on both kinds 
of networks. At low noise intensities the random networks are more robust against losing the coherency but upon 
P^ I increasing the noise, at a specific noise strength the synchronization among the population vanishes suddenly. 

fvj . In contrast, on scale-free networks the global synchronization disappears continuously at a much larger critical 

' noise intensity respect to the random networks. 
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I. INTRODUCTION 



(~| ■ Collective behavior in a population of individuals is of great importance in many areas of physics, biology, social sciences 

' and many other disciplines|l]J,|2,l3|]. One of the most celebrated collective behaviors is the case of phase synchronization among 
^ . a population of interacting self-sustained oscillators, in which all members tend to oscillate coherently with more or less the 
^ same phase. A population of coupled phase oscillators with mutual sine interaction between the pairs, is known as the Kuramoto 

1 K . model 1^. This model is introduced by Kuramoto Q and has extensively been investigated by many authors (see laSlMl and 
Ci • references therein). 

c/2 ] Considering a network of coupled rotors (phase oscillators), many factors such as couplirig strengthQlOll. time delayed inter- 

actions lUin . individual frequency distributionlalTtl, topology of network Ill2fl and noises||^[i3|l, affect the path toward the full 

2 [ synchronization. 
The synchronization of the deterministic Kuramoto model with random distribution of rotor frequencies and initial phases, has 

been studied on the scale-free networks by Moreno and Pacheco ll6\. There, it has been shown that the onset of synchronization 
occurs through a continuous transition at a small value of coupling with a critical exponent near 0.5. This resembles the mean- 
Q ' field behaviour, except that in the case of scale-free networks the critical coupling at which the rotors begin to get synchronized, is 
Q ' much smaller respect to all-to-all networks. They have also managed to find that in the complete synchronized state, dependence 
'~~' of the recovery time respect to the node degree is a power law function with exponent close to -1, which shows the robustness 

. [ of highly connected nodes (hubs) against perturbations. 
^ . The comparison between the synchronizability of the Kuramoto model on Erd6s-Renyi(ER) and scale-free networks has been 

fT^ ■ carried out recently by Gomez-Gardeii et al lUTlllSll . In these references, the authors have found that while the onset of global 
00 . synchronization occurs at smaller value of coupling for the scale-free networks, but tendency toward the global coherence grows 
CO ■ suddenly to larger extends for ER networks, at higher couplings. The reason for such behaviour is that the giant cluster in 
^N \ the heterogeneous networks (such as scale-free networks), originating from a central core of high connectivity (hub), grows 
■^ ' continuously by attaching the smaller clusters to it upon increasing the coupling constant. In contrast, for the homogeneous 
C^ networks, the evolution toward full synchronization would be boosted up by merging many small clusters spreading uniformly 
OO though out the networks, when the coupling is large enough. 

In this paper our focus is on the effect of white noise forcing on the synchronization of Kuramoto model for two types of 
networks: scale-free networks introduced by Barabasi and Albert(BA) Ill4ll and random networks [1 Si]. The paper is organized 
as follows: in section II, we introduce the stochastic Kuramoto model and express briefly some results for all-to-all networks, 
section III, is devoted to the simulation results and the conclusion will be presented in section IV. 
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II. STOCHASTIC KURAMOTO MODEL ON A COMPLEX NETWORK 

Consider a system composed of N rotors with the intrinsic frequencies, denoted by uji, on the top of a complex network 
consisting of N nodes. The stochastic Kuramoto model on such a network is described by the following set of equations: 

HO 

-^^^^ + >^J2 "'J ^"'(^J - Oi) + m(.t) , i = l---N, (1) 

j 

where 6i is the phase of the rotor on the node i, A is the coupling constant and aij is an element of the connectivity matrix, 
which takes the value a,., = 1 if z and j are linked together and otherwise, jy,; is the random force applying on the i-th rotor, 
and usually is chosen as a white noise with zero mean-value. The spatial-temporal correlation of such a noise is given by: 

Mt)Ti,it'))=2DS{t-t')S,j, (2) 

where D is variance of the noise. 

The stochastic Kuramoto model on an all-to-all network has been investigated analytically by Acebrdn et al ^, who have 
shown that taking a one peaked symmetrical frequency distribution /(w) = f{—'-o) for oscillators, there would be a critical 
coupling constant Ac = 2/[7r/(0)] above which the network begins to get synchronized. Near this critical point, the order 
parameter obeys a power law relation, namely 



-16(A-A.) ^3^ 



V 'rA4/"(0) ■ 

They have also shown that for a Lorentzian frequency distribution f{uj) = {j / tt) / {oj'^ + j'^), the incoherent solution is linearly 
stable for points (A, D) above the critical line D = —7 + A/2. In terms of coupling strength, this is also linearly stable for 
\< Xc = 2D + 27. 

In the next section we numerically integrate Eq.O on scale-free and random networks and compare the results. 

III. SIMULATION RESULTS FOR SCALE-FREE AND RANDOM NETWORKS 



To create a scale-free network with average connectivity (fc) = 2m, we use the BA algorithm lll4ll . In this procedure, starting 
from mo initial nodes all connected to each other, at each step one attaches a newly entering node to to < to-q elder ones such 
that the nodes with higher connectivity have larger probability (proportional to their degree) to get connected with this new one. 
Repeating this stages provides us with a network whose degree distribution obeys a power law function as P{k) ^ k '' with 
7 = 3. For producing ER network composed of A^ nodes and with the same average degree per node((fc) = 2m), it is enough to 
distribute Nm, edges between randomly chosen pair of nodes II15I1. 

In this work, we set m = 10 and select a delta function distribution for intrinsic frequencies, /(w) = S{uj — ujq). Changing 

the reference frame to a rotating one with rotation frequency ujq, enable us to set cjj = for all rotors in Eq.([T])- We also pick 

2 
77i(i) out of a box distribution in the interval ~g/2 < 77 < g/2, hence its variance is given by D = |j. Using Ito's formalism 

for integration of a stochastic function lll9ll . one obtains the following discrete equation from Eq.([T]i: 



i{t + dt) ^9,(t) + X 



j 



dt + j],{t)Vdt + 0{dr), (4) 



where in the Ito's picture, j]i (t) is evaluated at the initial point of the time interval [t,t + dt]. Time step, dt, is taken small enough 
to reduce the computational error The initial values of 9i are randomly drawn from a uniform distribution in the interval [— tt, it] . 
To characterize the global phase coherency, we define the following order parameter: 

1 ^ 
i=i 

which (• • •) means the averaging over different realizations of noise and initial conditions. In the stationary regime the time 
argument of r{t) can be omitted and one can replace the averaging over realizations by time averaging. The order parameter 
takes the value < r < 1, where r ~ corresponds to the disordered phase and r = 1 characterizes the full synchronized state. 



In the absence of noise for scale-free network, we found that the rotors get synchronized for a very small value of coupling 
around A = 0.03. We choose a large enough value for the coupling to make sure that the system is in the full synchronized state 
when the noise is turned off, then increase the noise intensity until the global coherency vanishes at the critical value of the noise 
strength, gc- 

In addition to the order parameter r, for better specifying the transition from coherency to decoherency, we introduce the 
Binder's forth cumulant which is defined as: 



" = i-^py2- (6) 

It is easy to see that in coherent phase, where (r) is nonzero, u takes the value 2/3 in the large N limit, while upon vanishing 
the global synchronization((r)=0) this quantity falls down to 1/3 in thermodynamic limit. The much smaller numerical errors 
in computation of the Binder's forth cumulant rather than the order parameter, makes it more advantageous for determination of 
the position and treatment of the coherency-decoherency phase transition. 

FigslT]and|2]represent the time dependence of order parameter for scale-free and random networks composed of A^ = lO** 
rotors, respectively. To derive these data, we put dt = 0.01, A = 0.2 and averagings have carried out over 100 different 
realizations of noise and initial phase configurations, for the noise intensities increasing from g — 2.0 by step /S.g ~ 2.0. From 
these figures, one finds that after about 500 time steps the system reaches the stationary for all noise intensities, and obviously 
the global coherency vanishes at larger coupling values for the scale-free network (around g = 10.0) respect to the random 
network (around g = 8.0). 

In what follows, to find the dependence of the order parameter as well as Binder's forth cumulant on the noise intensity, we 
fix the number of nodes to A^ = lO"', and the averagings are carried on lO"' time steps after skipping 2000 initial steps, where 
the system is surely in the stationary state. 

In Figs|3]and|4]we have depicted the order parameter, r, versus the noise intensity, g, for scale-free and random networks, for 
three coupling constants A = 0.1, 0.15, 0.2. Similar graphs for the Binder's forth cumulant, u, is represented in Figs|5]and|6] 

For comparison, the noise intensity variations of r and u have been depicted for both scale-free and random networks in 
Figs. 7a and 7b. 

By inspecting these figures one can extract two essential results: 
(i) Synchronizability of each kind of networks depends on the coupling strength, such that at small noise intensities the order 
parameter for random network falls more slowly than scale-free's, so it is more robust than SF network against the noise while 
at large noise intensities the situation is vice versa (see Fig. 7a). The critical noise intensity (gc) at which the transition from 
synchronized to un-synchronized state occurs is larger for the scale-free network. So the coherent state in SF network persists 
more against the noise than the random network with the same average degree and coupling constant. 

(ii) The coherency among the population of rotors destroys smoothly by increasing the noise intensity in the scale-free networks, 
while in the random networks the synchronization disappears by a sudden fall at the transition point. These behaviours are 
more apparent from the treatment of Binder's forth cumulant shown in Figs|5] |6]and 7b. Then the order-disorder transition in 
SF networks resembles the continuous transitions in equilibrium critical phenomena, while the transition in random networks is 
discontinuous-like. 

Referring to Gomez-Gardenes et a/ IUTlllSll . a nice explanation of our results are as follows. In homogeneous systems such as 
random networks, starting from the fully synchronized phase, when we turn on the noise, some incoherent clusters with more 
or less the same size begin to form. At low noise intensities, the size of these clusters are small and they are well separated, but 
by increasing the noise they get larger and connected to each other at intermediate noise intensities. At this point, the locally 
synchronized regions are not coherent anymore, so a big drop occurs for the order parameter This is much like the first order 
phase transitions in equilibrium statistical mechanics, where the ordered and disordered phases coexist at the transition point. 
On the other hand, the fully coherent state in the SF networks is founded around a core consists of nodes with high connectivity 
(hubs). When noise is applied on such state, the un-synchronized parts leave this giant cluster one by one, leading to continuous 
destruction of global coherency. 

IV. CONCLUSION 

In summary, we numerically investigated the stability of the global phase synchronized state in Kuramoto model on the top 
of scale-free and random networks, under white noise forcing on each oscillator. Our results emphasize on the fact that the 
stability of the synchronized phase is dependent on the noise strength, such that at low noise intensities the random networks are 
more stable against loosing the coherency, while at intermediate noise intensities, the coherency falls abruptly in such networks. 
However, in scale-free networks the coherency among the rotors decreases smoothly and also persists up to larger extends of 
noise intensity. Therefore, our findings confirm the picture presented by Gomez-Gardenes et a/ lll7lll8ll . that in heterogeneous 
networks the giant cluster formed around a core of hubs, grows (falls) continuously by increasing (decreasing) the coupling or 



by lowering (rising) the noise intensity. On the contrary, in homogeneous systems such as random networks, by increasing the 
noise intensity, the coalescence of un-synchronized clusters which are uniformly distributed over the network, results in a sudden 
fall in the global synchronization at the transition point. So the more complex is a system, more predictable it is. 

This work sheds more light on the different aspects of nonlinear dynamics on the top of homogeneous and heterogeneous 
network topologies and we hope that it promotes more researches on this very interesting problem. 
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FIG. 1 : Order parameter versus time for different noise intensities on tfie scale-free networlc. Tlie results are obtained for coupling constant 
A — 0.2 and the number of nodes is A^ = 10 . 
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FIG. 2: Order parameter versus time for different noise intensities on the random network. The results are obtained for coupling constant 
A — 0.2 and the number of nodes is A^ = 10*. 




FIG. 3: The order parameter versus the noise intensity for the scale-free network. The results are obtained for three coupling values A 
0.1, 0.15, 0.2 and iV = 10"* phase oscillators. 
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FIG. 4: The order parameter versus the noise intensity for the random network. The results are obtained for three coupling values A 
0.1, 0.15, 0.2 and iV = 10'' phase oscillators. 
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FIG. 5: The Binder's forth cumulant (Eq[6l( versus the noise intensity for the scale-free network. The results are obtained for three coupling 
values A = 0.1, 0.15, 0.2 and A'' = 10* phase oscillators. 
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FIG. 6: The Binder's forth cumulant (Eq|6l versus the noise intensity for the random network. The results are obtained for three coupling 
values A = 0.1, 0.15, 0.2 and N = lO'' phase oscillators. 
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FIG. 7: Noise intensity dependence of (a) Order parameter and (b) Binder's fortli cumulant for the scale free and random network at A = 0.1. 



